1 Introduction

The purpose of this document is to provide a reproducible record of all analyses and figures in the main article. Specifically, this document reproduce the same main analysis as in the main article, but excluding all extinct species from the analysis, as the inclusion of the many 0s, may change the results of the analysis.

2 Load datasets, Data wrangling and Imbalance calculation

3 Biomass

Let’s have a look at the biomass dynamics in the different environmental treatments.

3.0.1 tot biomass plot

Figure 1 : Community total biomass during the experiment in different environmental treatments. Different color represent richness levels.

4 Main Results

We now look at the main results of the experiment. We are going to look first at the effect of richness, temperature and nutrients on community temporal stability. Then, we are going to look at the relationship between divergence (original response diversity metric) and temporal stability. Finally, we are going to look at the relationship between imbalance and temporal stability.

In the whole analysis, we calculated the temporal stability of total community biomass as the inverse of the coefficient of variation (ICV) (i.e. \(\frac{\sigma}{\mu}\)).

4.0.1 Effect of T, N and R

Figure 2: Effects of richness (a), temperature (b), and nutrients (c) on community total biomass temporal stability.

4.0.2 Effect of Divergence

We look at the relationship between divergence (our original response diversity metric) and stability

Figure 3: Relationship between Divergence and temporal stability of total community biomass.

5 Effect RD

We are now going to look at how imbalance affected temporal stability of total community biomass. We are going to look at the relationship between fundamental imbalance (so based only on species response surfaces measured in monoculture), an realised imbalance (measured accounting for species contribution to balance).

This is fundamentally testing our most important hypothesis.

Figure 6: Effects of fundamental and realised imbalance on total community biomass temporal stability.

6 Linear models

6.1 Dependence of Imbalance on Temperature

Since imbalance depends on the slopes of species’ responses, and the slopes of species responses between two environmental conditions depend on the shape of species responses and the position of the optimum of the species response curve, we may expect that imbalance depends on temperature. For example, when the environment (i.e. temperature) fluctuates at higher mean temperatures, the species responses may be more similar (i.e. all species have the same sign of response), leading to less symmetrical response distributions, and thus higher imbalance.

We look now at how imbalance changes with temperature, in all possible combinations of species responses for each richness level and each environmental condition. We calculate the sum of slopes of all possible combinations of species responses, and then calculate the imbalance for each combination. We then plot the relationship between imbalance and temperature for each combination of species responses.

Figure 7: Relationship between imbalance and temperature for all possible combinations of species responses. The points represent the imbalance of each combination of species responses The points are colored by temperature regime, and the facets represent whether the combination of species responses was selected for the experiment.

We can see that imbalance increases with temperature, and that the relationship is the same in selected vs non selected community compositions. This suggests that imbalance is intrinsically dependent on temperatures, and that the species responses are more similar when the environment fluctuates at higher temperatures. Inevitably, this will affect the stability of the community, as we have seen in the previous analysis. Yet, it also means that the effect of imbalance on stability is not independent of temperature, or more generally, of the environmental conditions. Indeed, temperature and nutrients have been used to estimate species’ responses, and thus imbalance inherently depends on the environment.

Let’s look at the variance explained by imbalance and the environment on stability. Figure 8: Variance explained by imbalance and the environment on stability. The total variance explained by imbalance and the environment is shown in grey, and the unique variance explained by imbalance and the environment is shown in light grey.

95% of the variance explained by balance is shared with nutrients and temperature, which is not suprising considering how imbalance is calculated. However, there is additional variance explained by the environment that is not shared with balance. This suggests that the environment has an independent effect on stability, which is not mediated by imbalance. This may be related to the fact that imbalance was calculated using monoculture data collected at constant temperature, whereas in the community experiment, temperature was fluctuating, and species were interacting.

7 Model only looking at imbalance and richness, environmental conditions as random effects

We first investigate the effect of imbalance on stability while accounting for different environmental contexts. We fit a linear mixed-effects model including temperature, nutrient level, and composition as random factors to account for the different environmental contexts in which the imbalance was measured. This allows us to assess the effect of imbalance on stability while controlling for the influence of temperature and nutrient levels.

model_rev3 <-lmer(log10(stability) ~ log10(balance_f) + richness + 
       (1 + temperature + nutrients || composition),
     data = complete_aggr_2,
  REML = FALSE)
# check model's assumptions
check_model(model_rev3)
model check 1.

(#fig:model_check_int)model check 1.

Figure 9: Model check for the linear mixed-effects model with balance, richness as fixed effect. The model fits well the data.

Table 10: Linear mixed-effects model results for the effects of balance, richness on community stability. Estimates are presented with 95% confidence intervals and p-values.

Linear Regression Results
Predictor Estimate 95% CI t-value p-value
Intercept -0.16 -0.25, -0.08 -3.91 <0.001
log10(imbalance) -0.02 -0.04, 0.00 -1.58 0.12
richness



    2

    3 -0.06 -0.17, 0.06 -1.04 0.3
    4 0.04 -0.10, 0.17 0.58 0.6
Abbreviation: CI = Confidence Interval

Table 11: ANOVA table of the linear mixed-effects model with balance and richness as predictors of stability.

Term F Statistic DF p-value
log10(imbalance) 2.492960 1 0.114
richness 2.372282 2 0.305

7.1 Including the effect of environmental conditions on stability

To disentangle the shared variance between imbalance, temperature, and nutrients, we first removed the variance in temperature and nutrients that could be explained by imbalance. This can be done by regressing temperature and nutrients on imbalance and extracting the residuals, a method previously used to isolate independent effects among collinear predictors (e.g., Dormann et al. 2013; Graham 2003). The residual temperature and nutrient values thus represented variation independent of imbalance, allowing us to test their effects on stability without confounding influences of shared variance with imbalance. This approach enabled us to evaluate the relative importance of species richness and imbalance while accounting for environmental variability. By including residual temperature and nutrient values in the model, we ensured that our analysis focused on their independent contributions to stability.

7.2 Control effect of environment on balance

Now, we can build a linear model with imbalance, richness, the residuals of temperature and nutrients and their interaction as predictors of stability. We also use composition as random effect to account for some compositions being used in more than one treatment combination.

# check model's assumptions
check_model(model_no_random)
model check 1.

(#fig:model_check_int2)model check 1.

Figure 10: Model check for the linear mixed-effects model with balance, richness, and the residuals of temperature and nutrients as predictors of stability. The model fits well the data.

Table 12: Linear mixed-effects model results for the effects of balance, richness, and the residuals of temperature and nutrients on community stability. Estimates are presented with 95% confidence intervals and p-values.

Linear Regression Results
Predictor Estimate 95% CI t-value p-value
Intercept -0.17 -0.23, -0.10 -5.47 <0.001
log10(balance_f) -0.03 -0.05, -0.01 -3.16 0.002
richness



    richness3 - richness2 -0.02 -0.13, 0.09 -0.50 0.9
    richness4 - richness2 -0.04 -0.17, 0.09 -0.81 0.7
    richness4 - richness3 -0.02 -0.15, 0.11 -0.38 >0.9
resid_temp -0.04 -0.06, -0.02 -3.29 0.001
resid_nut 0.14 0.10, 0.19 6.27 <0.001
resid_temp * resid_nut -0.06 -0.13, 0.02 -1.46 0.15
Abbreviation: CI = Confidence Interval

The linear mixed-effects model was fitted by maximum likelihood (ML) with log-transformed stability as the response variable. The fixed effects included log10(balance_f), richness, the residuals of temperature (resid_temp), and nutrients (resid_nut), as well as their interaction term. The random effect was specified for the composition group, which accounted for potential variability between different composition categories.

Results mirror the previous analysis, showing that imbalance and the environmental variables have a significant effect on stability. The interaction between the residuals of temperature and nutrients is not significant when excliuding extinct species.

Table 13: ANOVA table of the linear mixed-effects model with balance, richness, and the residuals of temperature and nutrients as predictors of stability.

Term F Statistic DF p-value
(Intercept) 29.9437460 1 0.000
log10(balance_f) 9.9911933 1 0.002
richness 0.8136055 2 0.666
resid_temp 10.8050208 1 0.001
resid_nut 39.3231151 1 0.000
resid_temp:resid_nut 2.1207282 1 0.145

8 Asynchrony

Response diversity (one of the stabilisng effects captured by imbalance) has been suggested as a mechanism that promotes temporal stability of community biomass by promoting species asynchrony.

We thus calculated the asynchrony index suggested by Gross et al. 2014 to calculate the effect of asynchrony on temporal stability and to see how response diversity relate to asynchrony. The index ranges between -1 and 1, with -1 indicating perfect asynchrony and 1 being perfectly synchronous, and 0 indicating random variation.

8.0.1 Plot stability vs. Asynchrony Gross

Figure 11: Relationship between temporal stability and asynchrony (Gross) divided by nutrient level.

8.0.2 Plot Asynchrony Gross vs fundamental imbalance

Figure 12: Relationship between asynchrony (Gross) and fundamental imbalance divided by nutrient level.

8.1 Eveness

Evenness in species biomass has been identified as an important factor potentially influencing ecosystem stability Thibaut & Connolly 2013. In the context of our experiment, evenness in species biomass could help explaining why there is little difference between fundamental and realized imbalance. If evenness is high, then all species contribute similarly to total biomass. In this case, weighting for species-biomass contribution to total biomass (realized), should not fundamentally change the result, compared to an un-weighted (fundamental) measurement. Figure 13: Distribution of species evenness across experimental communities. The histogram represents the frequency of observed evenness values, while the red dashed line indicates the mean evenness (0.7). This highlights the central tendency of evenness across the dataset and its variation among communities.

Evenness was indeed generally high in our experimental communities, suggesting another potential factor reducing the potential difference between fundamental and realized balance.

9 Population stability

The relationship between community stability and the stability of the individual populations that make up the community is a key question in ecology. Importantly, ecosystem stability can result from low population stability, if populations fluctuate asynchronously, or from high population stability, if populations do not fluctuate much. Synthesis of the literature suggests diversity can have a positive or negantive effect on population stability Campbell et al 2010 and (Xu et al 2021)[https://onlinelibrary.wiley.com/doi/full/10.1111/ele.13777].

Theoretical work has suggested that community stability is a product of two quantities: the (a) synchrony of population fluctuations, and an average species-level population stability that is weighted by relative abundance Thibaut & Connolly 2013.

Critically, a imbalance value close to zero can result from high response diversity, but also from high population stability (population biomass does not change largely over time). We want to look now at whether our new metric of imbalance can capture these two stabilising mechanisms.

Thus, we calculate species-level population stability weighted by relative abundance and look at how it relates to ecosystem stability.

Figure 14: Relationship between log10 of population stability and log 10 of ecosystem stability.

Figure 15: Relationship between fundamental imbalance and population stability divided by nutrient level.

10 SEM

Finally, we use a structural equation model (SEM) to explore how stability is influenced by asynchrony, population stability, imbalance and, nutrient levels. In order to develop a hypothesis regarding the influence of stability, we have drawn on existing literature. This has enabled us to posit that stability is influenced by two key factors: asynchrony and population stability. In turn, these are influenced by balance and, in our particular case, by nutrient levels.

Path Coefficients (SEM, lavaan)
Response Predictor Estimate SE z-value p-value Standardized Estimate
stability asynchrony_Gross 0.136 0.023 5.833 0.000 0.348
stability Population Stability 0.449 0.053 8.444 0.000 0.573
asynchrony_Gross log10(imbalance) −0.093 0.034 −2.767 0.006 −0.185
asynchrony_Gross nutrients −0.211 0.023 −9.072 0.000 −0.503
Population Stability log10(imbalance) −0.073 0.011 −6.852 0.000 −0.292
Population Stability nutrients 0.136 0.008 17.854 0.000 0.653